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Electric current in semiconductors is considered as a function of time. Nonequilibrium fluctuations 
of the current are studied for strongly nonuniform samples. It is demonstrated that at special 
conditions the total electric current through a device can exhibit a negative fluctuation, when the 
qq ■ current is directed against the applied voltage. 
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The key to the understanding of the physical properties of semiconductor devices lies in the thorough analysis 
of equations modelling their behaviour [1]. The equations describing semiconductor transport characteristics are 
quite complicated, so that usually only stationary regimes can be satisfactory analysed. The situation becomes even 
more complicated when the distribution of charge carriers is essentially nonuniform. At the same time, strongly 
i nonequilibrium and nonuniform semiconductors can display peculiar features that cannot appear in stationary and 
^\ ' uniform samples [2-4]. 

In this communication we study nonequilibrium transport properties of semiconductor samples with a strongly 
nonuniform distribution of carrier densities. We show that at some special conditions there can appear a negative 
electric current, that is, the current directed against the applied voltage. The occurrence of this negative current is 
a transient effect that can be realized at the initial stage of a nonequilibrium process, while the carrier densities are 
yet essentially nonuniform. As soon as the charge distribution becomes more uniform, the current returns to normal. 
Thus, the negative-current effect is displayed as a kind of fluctuation happening because of strong nonuniformity of 
carriers in a nonequilibrium semiconductor. 

The charge carriers are characterized by the densities pi = pi(r,t), with i = 1,2, which are functions of the 

space variable r and of time t. Let p\ > and pi < 0. Carrier transport is commonly described in terms of the 
semiclassical drift-diffusion approximation which is the basis of semiconductor device models [1] . This approach rests 



q , on the continuity equations 



5-H ' 

and on the Maxwell equation 



with the electric-current density 



Vj< +HPi = (1) 



e VE=Mpi +P2), (2) 



j i~ PiPi E —Di V p tl (3) 



where 7.; is a relaxation width, e is the dielectric permittivity, p\ > and P2 < are the carrier mobilities, and Di is 
a diffusion coefficient. 

Adding to the sum of the drift-diffusion currents in eq.(3) the displacement current, one gets the total current 
density 

Jtot=3i + J2+^- (4) 
The total current across a device is given by the integral 
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J(t) = J 7 tot (r,t)dr, (5) 

with the integration over the considered sample. 

Consider a plane device of the width L and area A, being biased with an external voltage Vq. Then, instead of the 
vector r, we have one space variable x £ [0,L]. Define the effective transient time 

L 2 1 
T o = —FT, M = o (Mi + IM2I), 

where p is the average mobility. It is convenient to scale physical characteristics so that to deal with dimensionless 
quantities. To this end, we shall measure in what follows the space variable x in units of L; time t, in units of r ; and 
other physical quantities, in the corresponding units written below: 

Pa = -^jr, Qo = eAEq, E = -j-, 



7o = — , D = pV , 
to 
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For the plane geometry considered, eq.(l) reduces to 



dpi 
dt 



d 



d 2 Pl 



and eq.(2) reads 



dE 
dx 



dx 2 



4ir( Pl +p 2 ), 



where < x < 1 and t > 0. 

Assume that at initial time the sample is prepared with a nonuniform distribution of charge carriers 

Pi(x,0) = fi(x). 



(6) 



(7) 



(8) 



There exist several ways of forming samples with nonuniform carrier densities. This can be done, e.g., by means of 
special procedures of growing a semiconductor with doped layers incorporated during the growth [1] . Another way is 
to subject the sample to the influence of irradiation by particle or laser beams [3]. 

Let the considered device be biased with an externally applied voltage Vo, which, in the dimensionless units accepted, 
implies the condition 
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E(x,t)dx = 1. 



For the electric field we may derive from eq.(7), with condition (9), the functional 



Q(x,t) 



Q(x, t)dx 



E(x,t) = 1 + 4tt 

of the charge densities entering through the notation 

Q(x,t)= f [ Pl (x',t) + P2 (x',t)}dx'. 
Jo 

The total electric current (5) across the plane device can be written as 



(9) 



(10) 



J{t) = 

where condition (9) is taken into account and 



\ji(x,t) +j 2 (x,t)]dx, 



(11) 
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From eq. (11) we find 



J(t)= / [fiip-i(x,t) + H2p2{x,t)]E(x,t)dx+ 
Jo 

+ D 1 [ Pl (0, t) - P2 (l, t)] + D 2 [ P2 (0, t) - p 2 (l, t)] . (12) 

To understand the time behaviour of the electric current, we need to analyse expression (12), in which pi and E are 
defined by eqs.(6) and (7), with the initial condition (8) and the voltage condition (9). Such an analysis can be done 
resorting to computer calculations. However, before doing this, it is very useful to derive an approximate analytical 
solution clarifying the physics of processes. 

An approximate solution to the system of eqs.(6) and (7) can be obtained employing the method of scale separation 
[5-7] which is a generalization of the Krylov-Bogolubov averaging method [8] and the of guiding-centre approach [9]. 

The first step in the method of scale separation [5-7] is the classification of solutions onto fast and slow. In the 
present case, this is done as follows. Notice that the electric field is expressed as the functional (10) containing the 
carrier densities integrated over space. Thence, E should vary in space slower than pi. On the other hand, the voltage 
condition (9) shows that the electric field, being averaged over space, does not depend on time. This means that E can 
be regarded as a slow function in time. Therefore, the electric field can be classified as a slow solution, with respect 
to both space and time, as compared to the carrier densities. This allows us to consider eq.(6) for a fast solution pi 
treating there E as a space-time quasi-integral. Solving eq.(6), the found pi is to be substituted into functional (10), 
which results in an equation for E to be solved iteratively. 

When solving eq.(6), we continue pi outside the region < x < 1 by defining pi as zero for x < and x > 1. As a 
result we find 

/ + OC 
Gi{x-x',t)fi{x')dx', (13) 
-OO 
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where the Green function is 

ClAnr i\ = , , 

2y/^D~t \ ADit 

As an initial condition in eq.(8), it is reasonable to accept the physically realistic case of the Gaussian distribution 

fi{x) = — exp <j } (14) 

centered at the point a, € (0, 1) and in which 



Zi r I 2b t 
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J fi(x)dx, Z % = j exp <| - ^ ) dr. 



Then, from eq. (13) we obtain 

As is evident, at t = solution (15) reduces to the initial condition (14). And at t — > oo, we have 

lim pi(x,t) = 0, lim E(x,t) = 1. (16) 

t — >oo t — *oo 

Now, we have in hands enough information to try to find out whether a negative electric current could really appear, 
that is, an electric current directed against the applied voltage Vq. Take, for concreteness, that Vo > 0. Then we need 
to find when J(t) < 0. 

First thing, we immediately may notice, is that if Pi(x, t) is uniform in space then eq. (12) tells us that the electric 
current is a positively defined quantity J{t) = p.\p\ + \i2P2 > 0. So, if the electric current has a chance to become 
negative, this can happen only when pi(x,t) is nonuniform in space. 

From solution (15) it follows that even if the carrier densities are nonuniform at the initial time, nevertheless, they 
tend to become uniform with time, changing to zero, as t — > 00. Consequently, if a negative electric current could arise, 
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this might happen only at the initial stage when t< 1 and pi{x,t) are yet essentially nonuniform. Thus, a negative 
electric current, if any, can occur only as a principally transient effect at time t ^ I; and the space nonuniformity of 
the initial carrier density pi(x,0) — fi(x) is a necessary condition for this effect. 

Any initially nonuniform density, according to form (15), spreads with time because of diffusion and relaxation 
effects. However, these processes need some time to smooth pi. Hence, there always can be found such a time t <C 1 
when pi is yet nonuniform. But diffusion and relaxation shorten the time of a negative-current fluctuation, if such 
arises. Therefore, the conditions favoring the appearance of the noticeable negative-current effect are A < 1 and 
It < 1. 

After realizing what are the necessary and favoring conditions for the transient effect of a negative-current fluctu- 
ation, let us elucidate sufficient conditions for the inequality J(t) < 0. Since the latter can happen only at t <C 1 and 
with a nonuniform density of carriers, let us take t — and the maximally nonuniform initial density 

P i(x,0) = f i (x)=Q i 5(x-a i ), (17) 

resulting form eq.(14) under 6j — > 0. Then eq.(12) yields 

J(0) = mQtEiauO) + p 2 Q 2 E(a 2 , 0). (18) 

And for the electric field (10) we get 

E(x, 0) = 1 + 4ttQi [ai - 0(ai - x)} + AnQ 2 [a 2 - 6(a 2 - x)] , (19) 
where Q(x) is the unit-step function. From eqs.(18) and (19), we find that the inequality J(0) < acquires the form 

MiQi jl + 4ttQi ^ai - +4irQ 2 [a 2 - 6(a 2 - &i)]| + 

jl + 4^Q 2 (a 2 - + 47rQ! [ai - 0(ai - a 2 )]| < 0. (20) 

Inequality (20) is a sufficient condition for the occurrence of the transient effect of a negative electric current. 

There exist a number of cases when inequality (20) holds true. To prove that such cases do exist, consider a 
particular example when a\ = a 2 = a. Then, from eqs.(18) and (19), we derive that J(0) < means that 



• ' 1 1 • \-Q {<i ^ ) s >- (21) 

where 

Q = Qi+Q 2 . (22) 
Depending on whether the charge (22) is positive or negative, we have from (21) either 



or, respectively, 



«>\ + ^ «<0). (24) 



Remembering that < a < 1, we see that both eqs. (23) and (24) can be satisfied only when 

IQI > ^- (25) 

Eqs. (23), (24) and (25) are sufficient conditions for the appearance of a negative electric current at the initial stage 
of the process starting with the nonuniform distribution (17). 

To confirm the above analysis, we have solved eqs. (6) and (7) numerically and calculated the electric current (12) 
as a function of time. We accept the case favoring the appearance of the negative current, when Di -C 1 and 7$ <C 1. 
The initial condition (8) is taken in the Gaussian form (14). The voltage integral (9) plays the role of a boundary 
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condition for the electric field. For the carrier densities we may accept the Neumann or Dirichlct boundary conditions 
[1]. We considered both of them and found that the general behaviour of solutions is practically the same, but the 
calculational procedure is less stable in the case of Dirichlet conditions. Therefore, for the reason of stability, we 
preferred the Neumann boundary conditions. The numerical solution of eqs.(6) and (7) was accomplished using the 
third-order Rusanov finite-difference scheme [10]. 

The results of numerical calculations are in complete agreement with the necessary and sufficient conditions, ob- 
tained above, for the appearance of a negative electric current. Consider, e.g., the case of unipolar or majority-carrier 
devices such as metal semiconductor field-effect transistors, Schottky varactor diodes, and GaAs transferred electron 
devices [1] that are modelled by a single species set of the basic transport equations. In the case of one type of charge 
carriers, all processes are absolutely similar with respect to whether this is positive or negative charge. Let us take 
Qi = Q ^ and Q 2 = 0. Also, for brevity, we shall write a\ = a and b\ = b. The time behaviour of the electric 
current (11) is presented in Figs.l to 5. 

Fig.l demonstrates that if condition (25) does not hold then the negative electric current does not appear. In Fig. 2, 
the charge satisfies condition (25), and the negative current arises provided condition (23) is valid. Fig. 3 emphasizes 
the necessity of an essentially nonuniform distribution of the carrier density at initial time. Fig. 4 illustrates that 
a nonuniform initial density is necessary but not sufficient: to produce a negative current, the location of this 
nonuniformity is to satisfy condition (23). Fig. 5 stresses the importance of the latter condition showing that no 
negative current appears, even for a large nonuniform initial density, if this condition does not hold true. 

Now let us discuss the following two important questions: (i) How the studied electric current J(t) can be measured? 
(ii) What could be possible physical applications of the discovered effect? 

The quantity that can be measured directly, as is known [1], is the density of the total current across the ends of a 
device. That is, the directly measurable quatities are the density of current across the left end, 

J(0,t) = j tot (0,t), (26) 

and the current density across the right end, 

J(l,t) = j tot (l,t). (27) 

How the measurable end currents (26) and (27) are related to the studied total current (11)? 

According to the detailed analysis made above, the negative-current effect can arise as a transient phenomenon 
at the initial stage of the time evolution. Consider the relaxation process in the interval of time < t < t , with 
to <C 7j rl - In the case when ji -C 1, we may have t 3> 1, so that the considered time interval covers all stages of 
interest. For this time interval, the evolution equations, following from (1) and (2), in the plane geometry studied are 

dpi dji „ dE , , 

Basing on these equations, it is straightforward to derive that 

^3tot(x,t) = 0, (28) 

where j to t is the total current density (4). Eq.(28) tells us that j to t(x,t), actually, does not depend on the space 
coordinate x. If so, then it is obvious that 

J{t)= jtot(x,t)dx = jtot{x,t). (29) 
Jo 

Consequently, the end currents (26) and (27) are the same as the total current (29), 

J{0,t) = J(l,t) = J(t). (30) 

Thus, the studied total current J(t) is itself a directly measurable quantity. 

In numerical calculations, it is admissible to study any of the currents in eq.(30). However, the integral representa- 
tion (5) or (29) is more suitable for analytical investigation. This is why we dealt with this representation for deriving 
the explicit conditions under which the negative-current effect could arise. 

It is worth emphasizing that in all our consideration the total current density (4) always was taken in the complete 
form including the displacement current. The latter nowhere was neglected. The fact that, passing from eq.(5) to 
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eq.(ll), the displacement current formally disappeared in the integrand is due to the advantage of using the integral 
representation (5) for the total current and because of the voltage condition (9). Really, starting from eq.(5), we have 



J jtot{x,t)dx = J j{x,t)dx + ^- J ^E(x,t)dx, 



where 



j(x,t) = j 1 (x,t)+j 2 (x,t). 

Invoking the voltage condition (9), we get 

f 1 d d f 1 

J -E{x,t)dx = — J E{x,t)dx = 0. 

Consequently, we obtain the identity 

/ jtot(x,t)dx = / j(x,t)dx, 
Jo Jo 

Certainly, eqs.(5) and (11) would not be identical if the applied voltage would not be constant. If instead of (9) we 
would have the voltage condition 

/ E(x, t)dx = <p(t), 
Jo 

with a time-dependent voltage ip(t) in dimensionless units, then instead of (11) we would get 

f 1 d 
J(t)=J j(x,t)dx + ^-ip(t). 

The fact that under an oscillating applied voltage electric current oscillates as well can surprise nobody. This is why 
we concentrated on the case of a constant applied voltage. In the latter case, strong fluctuations of electric current 
are not expected. And the possibility that electric current turns against the applied voltage seems surprising. 

The discovered unusual effect of negative electric current is of interest by itself, irrespectively of its feasible appli- 
cations. Nobody can guarantee in advance what could be the most important use of a new effect. Nevertheless, we 
may suggest some possible physical applications. 

Since the appearance of the negative electric current is closely related to the properties of charge carriers in 
semiconductor, one can, by measuring this current, extract information about the characteristics of the carriers, for 
example, about their mobilities. 

Another possible application of the effect has to do with the case when the initial charge distribution is formed 
by irradiating a semiconductor sample by ions. In such a case, the charge distribution is a Gaussian centered at the 
distance, from the irradiated surface, of a mean free path of ions. The occurrence of the negative electric current 
essentially depends, according to eqs.(23) and (24), on the location of the initial charge distribution. Therefore, the 
negative-current effect can be used for measuring the mean free paths of ions in irradiating beams. 

One more application is related to the fact that semiconductor devices work very often under the influence of 
radiation. For instance, such devices can be a part of electronic systems serving near atomic stations or in cosmic 
space. The influence of radiation can lead to the development of nonuniformities spoiling the work of semiconductor 
devices. This, in turn, can result in misfunctioning of electronic systems, which can provoke dramatic consequences. 
Before this happens, it would be desirable to get a warning that nonuniformities due to irradiation have reached the 
dangerous level. Then one could prevent dramatic accidents. As far as a growing nonuniformity in a device can yield 
the occurrence of negative electric current, this effect could be used as a kind of warning signaling that the situation 
is getting dangerous. If inside an electronic system acting under the influence of irradiation one incorporates specially 
tuned semiconductor devices, based on the effect of negative current, then such devices can play the role of controllers. 

It is possible to speculate more suggesting other ways of using the effect. However, we would like to stress it again 
that a new physical phenomenon is of value as such. All experience teaches us that the most appropriate usage of a 
newly discovered effect practically never can be predicted in advance. 

In conclusion, we have shown that in noncquilibrium semiconductors a negative electric current can arise, being 
directed against the applied voltage. This effect is principally transient and can appear only at the initial stage of the 
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process, due to a nonuniform initial carrier density. The necessary and sufficient conditions for the occurrence of this 
effect are derived. The solution of the transport equations has been done both analytically and numerically. 
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Figure Captions 

Fig.l 

The total electric current as a function of time for the initial carrier density with Q = l/2n and b = 0.05, at two 
different locations: a = 0.25 (curve 1) and a = 0.5 (curve 2). All quantities here and in other figures are given in 
dimensionless units explained in the text. 



Fig.2 

The total electric current as a function of time for Q = 0.5 and b — 0.05, at different initial locations of carriers: 
a = 0.1 (curve 1), a — 0.25 (curve 2), a = 0.5 (curve 3) and a = 0.75 (curve 4). 



Fig.3 

The time dependence of the total electric current for Q = 0.5 and a = 0.25, at different widths of the initial 
nonuniform density of carriers: b — 0.05 (curve 1), b = 0.1 (curve 2) and b = 0.5 (curve 3). 



Fig.4 

The total electric current vs. time for Q = 3 and b = 0.05, at different initial carrier locations: a = 0.1 (curve 1), 
a = 0.25 (curve 2), a = 0.5 (curve 3) and a — 0.75 (curve 4). 



Fig.5 

The total electric current vs. time for Q — 3 and a — 0.5, at different widths of the initial carrier distribution: 
b = 0.05 (curve 1) and b = 0.5 (curve 2). 
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